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Abstract: We present the first complete calculation of the curvature perturbation 
generated during the hybrid inflation waterfall, caused by the coupling of the water- 
fall field to a gauge field A whose kinetic function f 2 depends on the inflaton field. 
We impose an upper bound on the field W = fA which ensures that it has a negli- 
gible effect before the waterfall. We confirm the claim of Soda and Yokoyama, that 
the perturbation 5W generates a statistically anisotropic spectrum and bispectrum, 
which could easily be observable. We also discover a new phenomenon, whereby the 
time-dependent 'varyon' field W causes the inflaton contribution to vary during the 
waterfall. The varyon mechanism might be implemented also with a scalar field and 
might not involve the waterfall. 
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1. Introduction 

During the waterfall of hybrid inflation, the perturbation of the waterfall field gener- 
ates a contribution to the curvature perturbation. But its spectrum is proportional 
to k 3 which almost certainly makes it negligible on cosmological scales [[jj. For the 
waterfall to generate a contribution with a nearly flat spectrum, its onset should 
be modulated by (i.e. depend upon the value of) some field that is different from 
both the inflaton and the waterfall field, and whose perturbation has a nearly flat 
spectrum. 

Such modulation was first considered in using a scalar field (and further 
explored in || |], ^ [J). Then the Soda and Yokoyama [0] used instead a U(l) gauge 
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field. 1 In this paper we give the first complete treatment of that case. The gauge 
field A has a kinetic function f 2 that depends on the inflaton field. We impose a 
condition ensuring that W = fA has a negligible effect before the waterfall, and 
take into account both the perturbation of / and the possible time-dependence of 
W. We confirm the claim of [0] that the perturbation <5W can generate a statistical 
anisotropic contribution to (, at a level which could easily be observed. We also 
find a new effect, which is that the time- dependence of W can cause a significant 
variation in the (statistically isotropic) inflaton contribution to (. We dub this new 
effect 'the varyon mechanism' and note that the varyon field might not be a gauge 
field and might not act during the waterfall. 

We will take for granted the main ideas of modern cosmology described for 
instance in ||, and use the notation and definitions of ]10|, |9|. The unperturbed 
universe has the line element 

ds 2 = -dt 2 + a 2 (t)5 ij dx i dx j . (1.1) 

In the perturbed universe, we can choose a slicing (fixed t) and threading (fixed x), 
and write for a given quantity g(x, t) = g(t) +5g(~x, t). A different slicing with a time 
displacement St(x, t) gives a different perturbation 5g. If g is rotationally invariant 
we have to first order 

5g(x, t) - 5g(x, t) = -g{t)6t{x, t). (1.2) 

We will invoke this 'gauge transformation' without comment. In most cases g is 
homogeneous on one of the slicings. 

We denote the Fourier component by Sg^t) where k is the coordinate wavenum- 
ber. Cosmological scales (probed directly by the CMB anisotropy and galaxy sur- 
veys) range from k = ko = (aH) Q to k ~ e 15 k Q , where H = a/a and (aH) is 
evaluated at the present epoch so that a /k is about the size of the observable uni- 
verse. A scale is 'outside the horizon' if k < aH. Inflation corresponds to tn < 1 
where e# = —H/H 2 . Cosmological scales leave the horizon during inflation and enter 
the horizon during the radiation-dominated era leading to Big Bang Nucleosynthesis 
(BBN). 



2. The curvature perturbation ( 
2.1 Definition and 5N formula 

To define ( one smoothes the metric on a super-horizon scale, and adopts the co- 
moving threading and the slicing of uniform energy density p. Then [12[ 

C(x, t) = 5[\n a(x, t)} = 5[ln (a(x, t)/a(t))} = 5N(x, t), (2.1) 
1 This is extended to the non-Abelian case in [0. 
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where a(x, t) is the locally defined scale factor (such that a comoving volume element 
is proportional to a 3 (x, t)). The number of e- folds of expansion iV(x, i, £*) starts from 
a slice at time £* on which a is unperturbed ('flat slice') and ends on a uniform p slice 
at time t. Since the expansion between two flat slices is uniform, 5N is independent 
of t*. 

By virtue of the smoothing, the energy conservation equation is valid locally: 

pW = -3^|^(p(t) + P(x,t)). (2.2) 

In consequence, ( = during an era when P(p) is a unique function. The success of 
the BBN calculation shows that P = p/3 to high accuracy just before cosmological 
scales start to enter the horizon. Then ( has a time-independent value £(x) that is 
strongly constrained by observation. Within observational errors it is gaussian and 
statistically isotropic. Its spectrum is nearly independent of k, with |L3| 



V c (k) ~ (5 x 1(T 5 ) 2 (2.3) 
n(k) - 1 = d\nV ( /d\nk =~ -0.032 ± 0.012. (2.4) 

(The result for n(k) assumes that it has negligible scale dependence. It also assumes 
a tensor fraction r <C 10 -1 , which will soon be tested by PLANCK flT4"| .) For the 
reduced bispectrum |L3f /nl, current observation give |/nl| ^ 100 and barring a 
detection PLANCK will give |/nl| ^ 10. For / NL to ever be observable we need 
I/nl| > I- 

We will work to first order in £, so that 

C(x,t) = J ff(t)(Jt /p (x,t), (2.5) 

where is the time displacement from the flat slice to the the uniform-p slice. 
A second-order calculation of £ is needed only to treat very small non-gaussianity 
corresponding to |/nl| ^ 1- 

We adopt the usual assumption, whereby iV(x, t, t*) is determined by the values 
of one or more fields </>j(x, t), evaluated during inflation at an epoch t*: 2 

iV(x,t) = iV(^(x),^(x),---,t). (2.6) 



Defining the perturbations Sep* on a flat slice, one writes [II], [15 



((x,t)=^iV 1 (t)#;(x) + ^iV ij (t)^(x)#;(x) + ..., (2.7) 



2 



2 To be more precise, JV will depend also on some of the masses and couplings in the action, and 
it may depend too on the values of any fields with negligible dependence on x and t that have not 
time to reach their vacuum expectation values. That does not affect any of the following. 
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where a subscript i denotes d/d(f>* evaluated at the unperturbed point of field space. 
The 4>i are usually taken to be scalar fields, but it has been proposed |7|, [10] that 
some or all of them may be components of a vector field. 

On each scale k, the field perturbations are generated from the vacuum fluctua- 
tion at horizon exit and are initially uncorrelated. Ignoring scales leaving the horizon 
after Eq. ( |2.7| ) defines a classical quantity (. In general it depends on t, settling 
down to the observed quantity C(x) by some time tf. Since C(x) is nearly gaussian, 
one assumes that Eq. (|2.7|) is dominated by one or more linear terms involving nearly 
gaussian scalar fields. With chosen as the epoch of horizon exit for a scale k this 
gives 

V c (k, t)~J2 Nf(U(k),t)V s<t>m (k, U(k)) + . . . , (2.8) 

where the terms exhibited correspond to scalar fields, and the dots indicate vector 
field contributions [TD|. Each contribution is positive. 



2.2 Slow-roll inflation 

Slow-roll inflation invokes Einstein gravity, and one or more scalar fields with the 
canonical kinetic term. The fields have practically gaussian perturbations, with 
^5</>i* — {H/27r) 2 at horizon exit. During single-field slow-roll inflation, only the 
inflaton has significant variation. Its unperturbed value 4>{t) satisfies 

3#0~-V'(0),, (2.9) 

where the potential V satisfies 

6 = l -Ml{V'/Vf^e H <^\ (2.10) 
H<1, ri = M$V"/V,, (2.11) 

giving p = 3MpH 2 ~ V. 

The perturbation 5(j)* generates a contribution Since 0(x, t) is the only time- 
dependent field, the effect on N of its perturbation 5(j)* (x) can be removed by the time 
shift St(x) which makes 0* homogeneous, which means that ^ is time- independent. 
At first order, 

C*(x) = -(tf/0)50»(x), (2.12) 

and 

7? (k)~—!—(—\ (2.13) 

C " V ; 2eM 2 \2n J K J 

n^k) - 1 = dV c Jdlnk = 2r l -Qe~2r l (2.14) 
where the right hand sides are evaluated at horizon exit. The second equality of 



Eq. ( 2.14 ) is appropriate for small- field models [16] and it applies to the standard 
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hybrid inflation which we are going to consider. The contribution of to |/nl| is 

Q<io- 2 . 

For multi-field slow-roll inflation, where two or more fields have significant varia- 
tion during inflation, Eqs. (|2.12|) - (j2.14 ) refer to the contribution of the field pointing 
along the inflaton trajectory at horizon exit. Field perturbations orthogonal to the 
(single- or multi-field) trajectory give no contribution to £ at horizon exit, but may 
contribute later. (This may occur during slow-roll inflation in a multi-field model, 
or during the waterfall, or after inflation through a curvaton-type mechanism.) We 
therefore have 

r u (k)<V c (k), (2.15) 

where V$(k) ~ (5 x 1CT 5 ) 2 is the observed quantity. The tensor fraction therefore 
satisfies 

r < 16e, (2.16) 

with e evaluated when fc leaves the horizon. 3 This leads []I6] to what has been called 
the Lyth bound, on the variation A0 of the inflaton field after k leaves the horizon, 
r < 10 _1 (A0 /Mp) 2 . For the tensor fraction to be detectable in the foreseeable future 
one needs r > 10 -3 [18|], which is impossible in a small-field model (A<p < 10 -1 M P ). 



3. The model 



3.1 Hybrid inflation 

We are interested only in the era starting with horizon exit for /c and ending with 
the onset of the waterfall. The relevant part of the action is taken to be 



S 



d X\ 



-fj 



-MiR 



1 



K, + AV(0) + ^m 2 ( 



m 2 ((f), A) = h 2 A 2 + g 2 (j) 2 - m 2 



A) X 2 + \\X* 



,(3.1) 

(3.2) 
(3.3) 



with F^ v 

we set / : 



: d^By — d v B^ and a U(l) gauge field. To fix the normalization of /, 
1 at a time t w just before the waterfall begins. 
Following [ID| we use the gauge with Bq = diBi = 0, and work with Ai = Bi/a 
which is the field defined with respect to the locally orthonormal basis (as opposed 
to Bi which is defined with respect to the coordinate basis). The raised component 
is A 1 = Ai (as opposed to B % = Bi/a 2 ). We also define the canonically normalized 
field W ee fA. The waterfall field x is supposed to be the radial part of a complex 
field which is charged under the U(l), generating the first term of Eq. (|3~ 



3 This follows from the definition r = Vh(k)/V^(k) with k ~ fc , and the prediction Vh(k) = 
(8/Mp)(iJ/2-7r) 2 with H evaluated at horizon exit. 
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We are going to impose Eq. ( ft. 12 ), which ensures that W has a negligible effect 



before the waterfall. Then, assuming suitable values for the parameters and field 
values, Eq. ( ft.l| ) gives what has been called |]J standard hybrid inflation |19|, At 



each location, the waterfall begins when m 2 (0, A) falls to zero. Before it begins, the 
waterfall field x vanishes up to a vacuum fluctuation which is set to zero, and we 
have slow-roll inflation with 

V = V + AV(<f>) ~ V Q . (3.4) 

We will take H to be constant which is typically a good approximation. In contrast 
with 0, we will not assume AV(4>) oc <fi 2 . 

During the waterfall, \ moves to it's vev and then inflation ends. We will assume 
that the duration of the waterfall is so short that it can be taken to occur on a 
practically unique slice of spacetime. This requires |]J m 3> H and H < 10 9 GeV. 

3.2 Field equations with / oc a a 

To work out the field equations, most previous authors have taken /(0(x, £)) to be 
a function only of time with / oc a a (t) (see however Taking spacetime to be 

unperturbed, the action (|3.1|) then gives for the unperturbed fields 



0(t) + 3Hj>(t) + V = (3.5) 

W(t) + 3HW(t) + /i 2 W(t) = 0, (3.6) 

where 

fi 2 = H 2 (2 + a)(l -a). (3.7) 

By virtue of the flatness conditions (|2.10|) and ( |2.11| ), the first expression is 



expected to give the slow-roll approximation (|2.9| ) more or less independently of the 
initial condition. Similarly, the second equation is expected to give the slow-roll 
approximation 3HW ~ — /i 2 W if \fi\ 2 H 2 . This condition is assumed because 
the analysis would otherwise become much more complicated. It is equivalent to 
a ~ 1 or —2. 

The first order perturbations satisfy 

64> k (t) + 3H6<j> k (t) + [{k/af + V"((f>(t)] 5<f> k (t) = (3.8) 
5W k (t) + 3F5W k (t) + {{k/a) 2 + /i 2 )5W k (t) = 0. (3.9) 



Keeping only super-horizon scales, Eqs. ( |3.6| ) and ( |3.9| ) and give 

3i7W(x, t) ~ -fx 2 deltaW{x, t). (3.10) 

The effect of the metric perturbation (back-reaction) on these equations vanishes 
in the limit where <fi and W are constant |9|, ^2[. The assumption of unperturbed 
spacetime is therefore expected to be a good approximation. 
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In terms of W, the coupling h 2 A 2 x 2 becomes h 2 W 2 x 2 , where h = h/f. We are 
setting / = 1 when the waterfall begins at t = t w . To generate <5W from the vacuum 
fluctuation, one assumes that W is a practically free field while cosmological scales 
leave the horizon, corresponding to h «C 1, or h <C e~ N ^ a where N(k) is the number 
of e-folds of inflation after horizon exit. With a ~ 1 this would make h too small to 
have a significant effect. One therefore assumes a ~ — 2. 

The simplest supersymmetric hybrid inflation model |23| has AV increasing log- 
arithmically. Then f(<f>) increases exponentially. The same behaviour holds for the 
non-hybrid model with the full potential V ~ 3MpH 2 oc <p 2 . It might be reasonable 
in string theory f24 |, and which could correspond to an attractor [p5 |. 



3.3 Field equations with /(0) 

In this paper we recognise that / is supposed to be a function of the inflaton field 
<ft, while retaining the assumption / oc a a for the unperturbed quantity. Using the 
slow-roll approximation with a = —2 we have 

tf/=^=|:^^=-2 (3.H) 

d(p da at d(p V2eM P 

Since / is a function of cf>, the term —\f 2 F IMU F llu in the action couples and W 
so that the right hand sides of Eqs. (|3.5| ), ( |3.6|) , (|3.8|) and (|3.9| ) are non-zero. We 
calculate them in the Appendix, and show that they are negligible if 

ep "2 ep "2 eU — 6eM| 

where pw is the energy density of W. We will assume this condition. It implies 
Pw ^ Pi which also ensures that W has a negligible effect during slow-roll inflation. 
From Eqs. ( |2.3| ) and ( |2.13| ) the condition corresponds to 

^<^(zm"\ (3 ,3, 



where k is the scale leaving the horizon at time t. 
3.4 The perturbation <5W 

The evolution equation for W(x, t) is the same as that of a free scalar field with 
mass-squared p 2 , and we are assuming \p\ 2 H 2 . Treating the Fourier component 
8W\f_(t) as an operator and assuming the vacuum state well before horizon exit, one 
finds well after horizon exit the approximately scale-independent vacuum expectation 
value 

3 2 2p2 
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where hats denote unit vectors. From Eq. ( 3.10 ), the operator 5Wk has almost 
constant phase which means that 5Wt can be treated as a classical quantity with 
this correlator. 

The decomposition 

W(x,<) = W(t) + <5W(x,t) (3.15) 

is made in some box of coordinate size L around the observable universe, with W(£) 
the average within the box. After smoothing on a cosmological scale k, the spatial 
average of (SW) 2 (evaluated within a region not many orders of magnitude bigger 
than the observable universe) is of order \n(kL)(H/2ir) 2 . We assume W(t) H, 
which is reasonable because W 2 (t) at a typical position is expected to be at least of 
order the mean square of (5W) 2 evaluated within a box with size M 3> L [^(J. 

Including both the inflaton and W and assuming that cubic and higher terms 
are negligible, Eq. (p.7[) becomes |L0| 



C(x,t) = C«^(x,t) + Cw(x,f) + l N 4>4>(t)[S<i>*(x)} 2 + N^M^Wt^pAd) 

i 

C*(x,t) = ^(t)<ty,(x), (3.17) 
Cw(x,t) ee £tf 4 (t)*W7(x) + i^^(t)(5W7(x)«5W7(x), (3.18) 

i ij 

where the subscripts on N denote partial derivatives evaluated on the unperturbed 
trajectory. We are assuming Eq. ( |3.12| ), which ensures that before the waterfall (w is 
negligible while £(x, t) is close to the time-independent quantity given by Eq. ( |2.12| ). 

4. Effect of the waterfall on £ 
4.1 End-of-inflation formula 

Let us denote the contribution generated during the waterfall by £ w . 4 To evaluate 
it, we assume that the waterfall happens very quickly so that it can be regarded as 
taking place on a single spacetime slice. Then [[1], Q 



C w (x) = H5t pp {y) = H 



^Pw(x) _ <W(x) 

. pit*) p(t+) 



H S -^^H6t pw . (4.1) 



In this equation, t pp (x) is the proper time elapsing between a uniform-p slice at time 
t w just before the waterfall and a uniform-p slice at time t + just after the waterfall, 
while t pw is the same thing with the final slice the waterfall slice itself. 

This end-of-inflation formula actually holds if the waterfall slice is replaced by 
any sufficiently brief transition from inflation to non-inflation. In |I[ it is invoked 



4 Notice that we are using the subscript w to indicate the waterfall era. The quantity £ w is the 
total change in £ during the waterfall, which as we shall see can be different from the contribution 
C w of <SW. 
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for the transition beginning during the waterfall, at the epoch when the evolution 
of x becomes non-linear. We are here applying it to the entire waterfall. It was 
first given |2| with A in Eq. ( |3.3| ) replaced by a scalar field. In [gj the slope of the 
potential in the A direction was assumed to be negligible corresponding to single- 
field hybrid inflation, and the same assumption was made in several later papers [0]. 
The assumption was relaxed in ||, §, corresponding to what has been called [|S| 
multi-brid inflation. Following |7[ we are here taking A to be the magnitude of a 



U(l) gauge field. One can also replace A by a non-Abelian gauge field [§, £7], || 

4.2 Waterfall contribution: general formula 

Instead of calculating ( w directly, we calculate 

C(x,t + )=Cw(x) + C0(x), (4.2) 

where C/>( x ) is given by Eq. ( |2.12|) . We do this first without specifying the function 
m 2 ((f), A) or the nature of A. We define (p w (A) by m 2 ((j) w ,A) = 0. (If this equation 
has more than one solution (p w (A), we choose one of them.) The waterfall occurs 
when 0(x, t) = 0„(x, t). 

If t/ P (x) is the displacement from the flat slice at t w to the uniform p slice at t + 
we have C(x, t+) = if<5t/ p (x). Making the good approximation Stf p = Stf w , where 
tf w is the displacement from the flat slice to the waterfall slice, we have 

0(x, t w + 8t fp (x)) = 4>{U) + S(f>(x, t w ) + </)(t w )5t/ p (x) (4.3) 
w (x, t w + 5t fp (x)) = <p w (t w ) + 5<p w (x, t w ) + <p w (t w )5t fp (x), (4.4) 

where 5(f) and 5<j) w are defined on the flat slice. For the unperturbed values this gives 
4>{t w ) = 4> w (t w ). For the perturbations it gives 

C(x, t+ ) = H5t fp (x) = g ^w(x,t w ) - ^(x,M (4 5) 

- <M*w) 

During hybrid inflation (f> < 0, and we need 4>(t w ) < 4> w (t w ), or the waterfall will 
never start. 

Now we invoke Eq. ( |3.3| ). Discounting the strong cancellation m 2 ~ h 2 A 2 it gives 

<Mx,t) = -(m 2 - h 2 A 2 (x,t)) 1/2 * ~ ~ i^M. (4.6) 
6/ g 2 mg 

In most previous work, A is taken to be a scalar field. For single-field hybrid 
inflation 0, |3|, W is supposed to be negligible. Then W has a practically time- 
independent value and the waterfall slice corresponds to simply 4>(x,t) = w (x). For 
two-brid inflation 0, [|, || , and A have equal status and the time-dependence of </>„ 
is significant. We have checked that in this case, Eqs. ( |4.5|) and (JO]) are equivalent 
to the result (4.1) given in j5]. 
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In our case A is the magnitude of a U(l) gauge field with the action Q3.1J) . Let 
us first follow by setting a = —2. From Eq. ( ft.lOj ), this makes W(^,t) time- 
independent. Then, if we ignore the perturbation 5f we have A(~x, t) oc 1/f oc a 2 . 
This gives Eq. ( |3.16|) for £(x, t + ), with 5 



where 



and 



C^(x,t + ) = ^, (4.7) 
Mx,t + ) = f^, (4.8) 



X = (4.9) 
M P mg^2e(t w ) 



Cw(x) = -(X/W 2 (U)) \W(t w ) ■ W(x,t w ) + ~5W(x,t w ) • 5W(x,t w )J . (4.10) 

(We display t w for future reference, even though we are for the moment taking W to 
be time- independent.) We need X < 1 for the waterfall to end. 

Ignoring the time- dependence of A we have Cvk( x , t+) = Cw( x ), which is the 
result obtained in 0. We see that Qw becomes bigger when the time-dependence is 
taken into account. The effect of the time-dependence of A was previously considered 
29] , who conclude that it decreases Cw by a factor e~ 2N ( k ) making it far too small 



in 



to have an observable effect. But the calculation of [2^] is not from first principles 
because A is treated scalar field. 

We have yet to include 5f. This has no effect on (w, but for it cancels the 
effect of A giving 

C*(x,t + ) = C*(x). (4.11) 

In the above we worked with (p w (A) defined by m 2 ((fi w ,A) = 0. That allows 
comparison with previous work where A is a scalar field, and with |7|, where A 
is the magnitude of a gauge field. But in the latter case the calculation becomes 
simpler if we use instead 4> W (W) defined by m 2 (0 w , (W j f {4>w)) = 0, because it is 
W that decouples from <p. With our current assumption /i 2 = 0, W vanishes which 
means that (f) w (W) = 0. Evaluating 5(f> w (W) we again arrive at Eqs. ( |4.8| ) and ( |4.11| ). 



According to this derivation, the factor 1 — X which was absent in [0] comes from 



5 We have checked that in this and the following cases, the third and fourth terms of Eq. (3.16) 
are negligible. 
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4.3 The varyon mechanism 

Allowing fi 2 ^ 0, we have W 7^ and hence <p w {W) ^ 0. This gives 6 



C(x,t H 



1 + 



6# 2 1 - X 



-i 



Cw(x) 

i -x 



(4.12) 



We see that C/>( x > is altered by the waterfall, which seems to contradict the state- 
ment in Section [2.2| , that the slow- roll inflation result C<^( x ) given by Eq. ( 2.12Q 
persists even after slow-roll ends. There is in fact no contradiction, because the 
presence of W(t) means that we are not dealing with exact slow-roll inflation. As 
a result, <f>(x.,t) is not the only time- dependent field, and the effect on iV(x, t, i*) of 
its perturbation 50* is not removed by the time shift which makes 0* homogeneous. 
Note that the perturbation 5W*(x) is in this context irrelevant. It is the quantity 
Wit) that is causing the effect. 

In our case W is the magnitude of a gauge field and its effect on ( occurs during 
the waterfall. But in general, any time-dependent field might do the same thing, i.e. 
have a negligible effect during slow-roll inflation but a significant one later owing to 
the time-dependence of its unperturbed part. The only exception is if the field is a 
slowly rolling scalar field; in that case its effect is just to slightly alter the direction 
in field space of the trajectory, which as we discussed in Section [2.2| will still give the 
(almost unchanged) the time-independent perturbation C/>( x )- We propose to call a 
field causing the new effect a varyon. 

The varyon mechanism can remove the bound ( j2.15|) and hence the bound ( f2.16| ) 
on the tensor fraction. That might allow an observable tensor fraction within a small- 
field inflation model. The possibility of avoiding Eq. fl2.16| ) was mooted in [[H], [32 



but they did not find a mechanism. One can easily avoid Eq. ( [2.16|) by abandoning 
the canonical kinetic term for the inflaton [|33] and we are now, for the first time, 
pointing to a possible mechanism with the canonical kinetic term. 

In our case, the varyon W has only a small effect except in the the very fine- 
tuned regime (1 — X) <C \fi 2 \/H 2 . Even if its effect is significant, and reduces 
(positive /i 2 ) so as to give r > 12e, it cannot make r big enough to observe because 
the end-of-inflation formula requires H < 10 9 GeV corresponding to r < 10~ 10 . It 
remains to be seen if a different varyon, perhaps a scalar field, can give a more 
interesting result. 

4.4 Anisotropic spectrum and bispectrum 

Since we are taking W 3> H, the linear term of (w dominates, leading to a spectrum 
of the form [0, |K| 



V c (k) = V c (k) 



I - .' ( A-k) 2 



(4.13) 



6 In pOl ) this formula is given with 1 — X incorrectly replaced by 1 + X. That makes the first 
term must be close to I whatever the value of X . 
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On cosmological scales, observation requires \/3\ < 10 1 , and barring a detection 
PLANCK will give \/3\ < 1(T 2 @. We therefore have P c (k) ~ V c (k) with (Eq. (p|)) 



7> c (fc) ~ (5 x lfr 5 ) 2 



With the parameters constrained to give /3 1, we have [0, |1C 

li 4 W 2 (iw) e(t fc ) 







m 2 g 2 e(t w ) 
^( fc ) n ^-2 

r „ 2 y i- 2 

1 + , (4.14) 



6# 2 1 - X 



where is the epoch of horizon exit for the scale k, and V^, is given by Eq. ( |2.13|) . 
The first line is the result of JTJ, who assumed V^(k) = V^(k). The first factor of 
the second line drops that assumption while the second factor takes account of the 
inhomogeneity of /(</>). The third line allows /i ^ 0. 

Including the quadratic term of ( w we reproduce the result of |7|] for the reduced 
bispectrum: 

/NL = /^(l + r ni (ki,k 2 ,k 3 )) (4.15) 

where 

/»> = -(A-^-(A.k^4-(k faxA.k.KA.M + 2 

ft 2 

/nl = 3 Y" ( 4 - 17 ) 
5. Conclusion 

Although there is so far no evidence for statistical anisotropy of the primordial cur- 
vature perturbation (, mechanisms have been proposed for generating it. Most of 
them invoke a vector field. 

One mechanism takes the vector field to be homogeneous during inflation, but 
causes significant anisotropy in the expansion (for a recent review of this ap- 
proach see ||36|| ). Then the perturbations of scalar fields generated from the vacuum 
fluctuation will be statistically anisotropic, and so too will be ( on the usual assump- 
tion that it originates from one or more of these perturbations. 

A different mechanism takes the inflationary expansion to be practically isotropic, 
but generates a vector field perturbation from the vacuum fluctuation f7|, [Hj (for the 
most recent paper on this approach see f3Tfl). 7 

In this paper we have given the first complete treatment of the version of the 
second mechanism proposed in |7|], which couples the waterfall field to a gauge field 



7 The use of a vector field to generate a contribution to C was first mooted in 
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A whose kinetic function f 2 depends on the inflaton. We have confirmed their claim 
that the statistical anisotropy could easily be big enough to observe, and we have 
also discovered a completely new effect; if W = fA is time-dependent it causes the 
usual inflaton contribution ^ to vary during the waterfall. This 'varyon' effect might 
still occur if W is replaced by a time-dependent (but not slowly rolling) scalar field 
and it might have nothing to do with the waterfall. 
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A. Equations of Motion for </>(x, t) and W(x, i) 



Extremizing the action in Eq. ( |3.1|) with respect to fields 0, and their derivatives 
we obtain field equations 

[d, + d, In y/=g\ d»<j> + V+ -ff'F^F^ = 0; (A.l) 
fa + d. In y/=j\ fF*" = 0, (A.2) 

where g = det(g^ u ) and /' = df /d(j>. Choosing the temporal gauge B = and a 



line element of the unperturbed universe in Eq. Ql.lD , one finds equations of motion 
for the fields 0(x, t) and B(x, t) 



+ 3H<f) - a- 2 V 2 + V = --ff'F^, 



H + 2-1 Bi- a - 2 V 2 B, ; = a- 2 2^j-d,B 



f 



(A.3) 
(A.4) 



Recasting the above equations in terms of W = /B / a and dropping gradient terms, 
one arrives at equations of motion for homogeneous fields <p(t) and W(t) 



+ 3H<j) + V 



f 



W+ \ H- j \ W 



(A.5) 



W + 3HW + ( 2H 2 - Hj - j ) W = 0, 



(A.6) 



where we also used H ~ 0. 
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Decomposing the field W(x, t) as in Eq. (|3.15|) and similarly the field 0(x, t), 
we find equations of motion for <50(x, t) and 5W(x, it) from Eqs. ( |A.3| ) and (|A.4j) . 
Keeping only the first order terms and switching to the Fourier space they become 



Jk + 3H5<p k + ( — + V" 



7 



W+ | H- L I w 



5Wl 



(A-7) 



5W k + 3H5W k +\2H 2 -H--- \ 5W k + ^5W k 




,2 ' 



w. 



(A.8) 



With our choice / oc a a for the unperturbed /, the above expressions become 



i k + 3#50 k + ( ^ + V" 1 50 k 



-2a 



'2eM P 
5W k + 3#5W k 
-aW 



W + i?(l-a)W 

- + ii 2 pw k 



5W k + H{l-a) <5W k - aW^50 k (A-9) 



2eM P 



A- 2 



i k + ^(l + 2«) 50 k + — 50 k 



(A.10) 



The energy density of the vector field in Eq. (|3.1|) is given by [[39] p#(x, i) = 
—pF^F^/A. From this it is easy to see that the background value of p_e(x, t) is 
given by 



| JET- ^ | W 



-# 2 I^ 2 . 
2 



(ah; 



The right hand side of Eq.( [A.5[ ) is negligible if pb satisfies Eq. ( |3.12[ ). We now show 
that the same is true of the right hand sides of Eqs. ( |A.9| ) and ( |A.10| ). At the epoch 
k ~ aH, the terms on the left hand sides are of order if 3 and Eq. ( |3.12| ) ensures that 
the right hand sides are indeed much smaller. At the epoch aH/k = exp(Nk(t)) ^> 1, 
the first term of each left hand side is negligible. The other two terms are of order 
\r]\ = \ V"\/3H 2 for Eq. Q and of order \vw\ = \^\/^H 2 for Eq. flPDP . Eq. ( gT| ) 
ensures that the right hand side of Eq. (K$) is negligible, and it ensures that the right 
hand side of Eq. ([A.10Q is negligible if also \rjw\ 10~ 5 . But the latter condition is 
irrelevant, because its violation makes the time-dependence of W (coming then from 
the right hand side) negligible. 
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